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Abstract. The quantification of synapses is instrumental to measure 
the evolution of synaptic densities of neurons under the effect of some 
physiological conditions, neuronal diseases or even drug treatments. 
However, the manual quantification of synapses is a tedious, error-prone, 
time-consuming and subjective task; therefore, tools that might auto¬ 
mate this process are desirable. In this paper, we present SynapGountJ, 
an Imaged plugin, that can measure synaptic density of individual neu¬ 
rons obtained by immunofluorescence techniques, and also can be ap¬ 
plied for batch processing of neurons that have been obtained in the 
same experiment or using the same setting. The procedure to quantify 
synapses implemented in SynapGountJ is based on the colocalization of 
three images of the same neuron (the neuron marked with two antibody 
markers and the structure of the neuron) and is inspired by methods 
coming from Gomputational Algebraic Topology. SynapGountJ provides 
a procedure to semi-automatically quantify the number of synapses of 
neuron cultures; as a result, the time required for such an analysis is 
greatly reduced. 


1. Introduction 

Synapses are the points of connection between neurons, and they are 
dynamic structures subject to a continuous process of formation and elim¬ 
ination. Pathological conditions, such as the Alzheimer disease, have been 
related to synapse loss associated with memory impairments. Hence, the 
possibility of changing the number of synapses may be an important asset 
to treat neurological diseases HZ]. To this aim, it is necessary to deter¬ 
mine the evolution of synaptic densities of neurons under the effect of some 
physiological conditions, neuronal diseases or even drug treatments. 

The procedure to quantify synaptic density of a neuron is usually based 
on the colocalization between the signals generated by two antibodies [1]. 
Namely, neuron cultures are permeabilized and treated with two different 
primary markers (for instance, bassoon and synapsin). These antibodies 
recognize specifically two presynaptic structures. Then, it is necessary a 
secondary antibody couple attached to different ffuorochromes (for instance 
red and green; note, that several other combinations of color are possible) 
making these two synaptic proteins visible under the ffuorescence micro¬ 
scope. The two markers are photographed in two gray-scale images; that, in 
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turn, are overlapped using respectively the red and green channels. In the 
resultant image, the yellow points (colocalization of the code channels) are 
the candidates to be the synapses. 

The final step in the above procedure is the selection of the yellow points 
that are localized either on the dendrites of the neuron or adjacent to them. 
Tools like MetaMorph [5] or Imaged m a Java platform for image pro¬ 
cessing that can be easily extended by means of plugins — can be used to 
manually count the number of synapses; however, such a manual quantifica¬ 
tion is a tedious, time-consuming, error-prone, and subjective task; hence, 
tools that might automate this process are desirable. In this paper, we 
present Synap Count an Imaged plugin, that semi-automatically quantifies 
synapses and synaptic densities in neuron cultures. 

2. Methodology 

SynapCountJ supports two execution modes: individual treatment of a 
neuron and batch processing. 

2.1. Individual treatment of a neuron. The input of SynapCountJ in 
this execution mode are two images of a neuron marked with two antibodies 
(an image per antibody), see Figure SynapCountJ is able to read tiff 
(a standard format for biological images) and lif files (obtained from Leica 
confocal microscopes) — the latter requires the Bio-Formats plugin [8]. The 
following steps are applied to quantify the number of synapses in the given 
images. 

In the first step, from one of the two images, the region of interest (i.e. 
the dendrites where the quantification of synapses will be performed) is 
specified using NeuronJ m — an Imaged plugin for tracing elongated image 
structures. In this way, the background of the image is removed. The result 
is a file containing the traces of each dendrite of the image. 

Subsequently, the user can decide whether she wants to perform a global 
analysis of the whole neuron, or a local analysis focused on each dendrite 
of the neuron. In both cases, SynapCountJ requires additional information 
such as the scale and the mean thickness (that is determined by the size 
of the subjacent dendrite) of the region to analyze (see Figure — these 
parameters determine the area of the dendrite avoiding the background (i.e. 
all the non-synaptic marking). 

Taking into account the settings provided by the user, SynapCountJ over¬ 
laps the two original images of the neuron and the structure of the neuron 
previously defined. From the resultant image, SynapCountJ identifies the 
almost white points (the result of green, red, and blue combination) as 
synaptic candidates, and it allows the user to modify the values of the red 
and green channels in order to modify the detection threshold (see Figure]^. 

Once that the detection threshold has been fixed, the counting process 
is started. Such a process is inspired by techniques coming Computational 
Algebraic Topology. In spite of being an abstract mathematical subject. 
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Figure 1. Neuron with two antibody markers and its 
structure. Left. Neuron marked with the bassoon antibody 
marker. Center. Neuron marked with the synapsin antibody 
marker. Right. Structure of the neuron. 



Figure 2. SynapCountJ window to configure the analysis. 

Algebraic Topology has been successfully applied in digital image analy¬ 
sis daiaE]. In our particular case, the white areas are segmented from 
the overlapped image, and the colors of the resultant image are inverted — 
obtaining as a result a black-and-white image where the synapses are the 
black areas. From such an image, the problem of quantifying the number 
of synapses is reduced to compute the homology group in dimension 0 of 
the image; this corresponds to the computation of the number of connected 
components of the image. 

Finally, SynapCountJ returns a table with the obtained data (length of 
dendrites both in pixels and micras, number of synapses, and density of 
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Figure 3. SynapCountJ window to modify the 
threshold of the red and green channels. Left. Win¬ 
dow to fix the threshold of the image. Right. Fragment of the 
neuron image with the synapses indicated as the red areas 
on the structure of the neuron marked in blue. Moving the 
scrollbars of left window, the marked areas of the image are 
changed. 


synapses per 100 micron) and two images showing, respectively, the analyzed 
region and the marked synapses (see Figure]^. 

2.2. Batch processing. Images obtained from the same biological experi¬ 
ment usually have similar settings; hence, their processing in SynapCountJ 
will use the same configuration parameters. In order to deal with this sit¬ 
uation, SynapCountJ can be applied for batch processing of several images 
using a configuration file generated from the analysis of an individual neu¬ 
ron. 

For batch processing, SynapCountJ reads tiff files organized in folders or 
a lif file (the kind of files produced by Leica confocal microscopes), and using 
the configuration file processes the different images. As a result, a table with 
the information related to each neuron from the batch (the table includes 
an analysis for both the whole neuron and from each of its dendrites) is 
obtained. 


3. Experimental Results 

The original aim of SynapCountJ was the automatic analysis of synaptic 
density on neurons treated with SB 415286 — an organic inhibitor of GSK3, 
a kinase which inhibition was proposed as a therapy in AD treatment |1] — 
such a treatment, as it was previously demonstrated, promotes synaptoge- 
nesis and spinogenesis in primary cultures of rodent hippocampal neurons 
and in Drosophyla neurons Eli. In this setting, a comparative study has 
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1 

Tracing N1: 
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71 
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Tracing N2: 
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35 
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Tracing N3: 
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Tracing N4: 
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25 
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116 

164 

6 

Tracing N6: 

468.8438 

23.4422 

26 
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38 
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116 
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69 
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116 

164 

10 

Tracing N10: 
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45 

73.4851 

116 

164 

11 

Tracing Nil: 

355.7054 

17.7853 

26 

146.1884 

116 

164 

12 

Tracing N12: 

905.3750 

45.2688 

45 

99.4063 

116 

164 

13 

Totai Neuron 

10474.9103 

523.7455 

479 

91.4566 

116 

164 



Figure 4. Results provided by SynapCountJ. Top. Ta¬ 
ble with the results obtained by SynapCountJ. Bottom Left. 
Image with the analyzed region of the neuron. Bottom Right. 
Image with the counted synapses indicated by means of blue 
crosses. 


been performed in order to evaluate the results that can be obtained with 
SynapCountJ. 

Primary hippocampal cultures were obtained from PO rat pups (Sprague- 
Dawley, strain, Harlan Laboratories Models SL, France). Animals were 
anesthetized by hypothermia in paper-lined towel over crushed-ice surface 
during 2-4 minutes and euthanized by decapitation. Animals were han¬ 
dled and maintained in accordance with the Council Directive guidelines 
2010/63EU of the European Parliament. Primary cultures of hippocampus 
neurons were prepared as previously described in m- Briefly, glass cov- 
erslips (12 mm in diameter) were coated with poly-L-lysine and laminin, 
100 and 4 /ig/ml respectively. Neurons at a 10 x 104 neurons/cm^ density 
were seeded and grown in Neurobasal (Invitrogen, USA) culture medium 
supplemented with glutamine 0.5 mM, 50 mg/ml penicillin, 50 units/ml 
streptomycin, 4% FBS and 4% B27 (Invitrogen, CA, USA), as described 
before in [I]. At days 4, 7 and 14 in culture a 20% of culture medium was 
replace by fresh medium. Cytosine-D-arabinofuranoside (4 /iM) was added 
to prevent overgrowth of glial cells (day 4). 

Synaptic density on hippocampal cultures was identifled as previously 
described in [I]. In short, cultures were rinsed in phosphate buffer saline 
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Figure 5. Quantification of synapses. Left. Manual 
quantification of synapses. Right. Quantification of synapses 
using Synap CountJ. 


(PBS) and fixed for 30 min in 4% paraformaldehyde-PBS. Coverslips were 
incubated overnight in blocking solution with the following antibodies: anti- 
Bassoon monoclonal mouse antibody (ref. VAM-PS003, Stress Gen, USA) 
and rabbit polyclonal sera against Synapsin (ref. 2312, Cell Signaling, USA). 
Samples were incubated with a fluorescence-conjugated secondary antibody 
in PBS for 30 min. After that, coverslips were washed three times in PBS 
and mounted using Mowiol (all secondary antibodies from Molecular Probes- 
Invitrogen, USA). Stack images (pixel size 90 nm with 0.5 fim Z step) were 
obtained with a Leica SP5 Confocal microscope (40x lens, 1. 3 NA). Percent¬ 
age of synaptic change is the average of different cultures under the same 
experimental conditions. As a control, we used sister untreated cultures 
growing in the same 24 well multi plate. 

A total of 13 individual images from three independent cultures has been 
analyzed. In Figurej^we can observe that using a manual method to identify 
and count synapses, we obtain a mean of 24.12 synapses in control cultures 
and 16.74 in treated cultures. The results obtained with SynapCountJ are 
similar, there is a mean of 26.03 synapses in control cultures and 16.50 in 
the ones which have been treated. 

Notwithstanding the differences in the quantification, in both procedures 
we obtain almost the same inhibition percentage, a 30.51% manually and 
36.61% automatically. This shows the suitability of SynapCountJ to count 
synapses, meaning a considerably reduction of the time employed in the 
manual process. Namely, the manual analysis of an image takes approxi¬ 
mately 5 minutes; of a batch, 1 hour; and, of a complete study, 4 hours. 
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Software Language Underlying Types of Technique for 

Technology Images detection 


Green and Red Puncta 

Java 

ImageJ 

tiff 

Colocalization 

Puncta Analyzer 

Java 

Image J2 

tiff 

Colocalization 

SynapCountJ 

Java 

ImageJ 

tiff and lif 

Colocalization 

SynD 

Mat lab 

Mat lab 

tiff and Ism 

Brightness 

SynPAnal 

Java 


tiff 

Brightness 


Table 1. General features of the analyzed software. 


Using SynapCountJ, the time to analyze an image is 30 seconds; a batch, 2 
minutes; and, a complete study, 6 minutes. 


4. Discussion 

Up to the best of our knowledge, 4 tools have been developed to quantify 
synapses and measure synaptic density: Green and Red Puncta m, Puncta 
Analyzer [19], SynD [13] and SynPAnal [3] — a summary of the general 
features of these tools can be seen in Table [H The rest of this section 
is devoted to compare SynapCountJ with these tools, the comparison is 
summarized in Table [2l 

There are two approaches to locate synapses in an RGB image either 
based on colocalization or brightness. In the former, synapses are identified 
as the colocalization of bright points in the red and green channels — this 
is the approach followed by Green and Red Puncta, Puncta Analyzer and 
SynapCountJ — in the latter, synapses are the bright points of a region 
of an image — the approach employed in SynD and SynPAnal. In both 
approaches, it is necessary a threshold that can be manually adjusted to 
increase (or decrease) the number of detected synapses; such a functionality 
is supported by all the tools. 

In the quantification of synapses from RGB images, it is instrumental to 
determine the region of interest (i.e. the dendrites of the neurons where the 
synapses are located); otherwise, the analysis will not be precise due to noise 
coming from irrelevant regions or the background of the image — this hap¬ 
pens in the Green and Red Puncta tool since it considers the whole image for 
the analysis. Puncta Analyzer allows the user to fix a rectangle containing 
the dendrites of the neuron, but this is not completely precise since some 
regions of the rectangle might contain points considered as synapses that do 
not belong to the structure of the neuron. SynD is the only software that 
automatically detects the dendrites of a neuron; however, it can only be 
applied to neurons with a cell-fill marker, and does not support the analysis 
from specific regions, such as soma or distal dendrites. SynapCountJ and 
SynPAnal provide the functionality to manually draw the dendrites of the 
image; allowing the user to designate the specific areas where quantification 
is restricted. 
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Software 

Detection of 

Threshold 

Batch 

Dendrites 

Density 

Export 

Save 


dendrites 


Process¬ 

length 







ing 






Green and 

Red Puncta 

Not used 

/ 






Puncta Ana¬ 
lyzer 

Manual ROI 

/ 




/ 


SynapCountJ 

Manual 

/ 

/ 

/ 

/ 

/ 

/ 

SynD 

Automatic 

/ 

/ 

/ 


/ 

/ 

SynPAnal 

Manual 

/ 


/ 


/ 

/ 


Table 2. Features to quantify synapses and synaptic 
density of the analyzed software. 


The main output produced by all the available tools is the number of 
synapses of a given image; additionally, SynapCountJ, SynD and SynPAnal 
provides the length of the dendrites; and, SynapCountJ is the only tool that 
outputs the synaptic density per micron. All the tools but Green and Red 
Punctua can export the results to an external file for storage and further 
processing. 

images obtained from 


Finally, as we have explained in Subsection 2.2 


the same biological experiment usually have similar settings; hence, batch 
processing might be useful. This functionality is featured by SynapCountJ 
and SynD, and requires a previous step of saving the configuration of an 
individual analysis. SynPAnal does not support batch processing, but the 
configuration of an individual analysis can be saved to be later applied in 
other individual analysis. 

As a summary, SynapCountJ is more complete than the rest of available 
programs. It can use different types of synaptic markers and can process 
batch images. Furthermore, a differential feature of SynapCountJ is that 
it is based on a topological algorithm (namely, computing the number of 
connected components in a combinatorial structure), allowing us to validate 
the correctness of our approach by means of formal methods in software 
engineering 


5. Conclusions and Further Work 

SynapCountJ is an ImageJ plugin that provides a semi-automatic proce¬ 
dure to quantify synapses and measure synaptic density from immunofluo¬ 
rescence images obtained from neuron cultures. This plugin has been tested 
not only with neurons in development, but also with the neuromuscular 
union of Drosophila; therefore, it can be applied to the study of images 
that contain two synaptic markers and a determined structure. The re¬ 
sults obtained with SynapCountJ are consistent with the results obtained 
manually; and SynapCountJ dramatically reduces the time required for the 
quantification of synapses. 
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As further work, it remains the tasks of improving the usability of the 
plugin and including post-processing tools to manually edit the obtained 
results. Additionally, and since the final aim of our project is the complete 
automation of the whole process, it is necessary a procedure to automatically 
detect the neuron morphology. 

Availability and Software Requirements 

SynapCountJ is an Imaged plugin that can be downloaded, together 
with its documentation, from http://imagejdocu.tudor.lu/doku.php? 
id=plugin:utilities: synapsescountj : start, SynapCountJ is open source 
and available for use under the GNU General Public License. This plugin 
runs within both Imaged and Fiji m and has been tested on Windows, 
Macintosh and Linux machines. 
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